Effective potential for quantum scalar fields in a de Sitter geometry 
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We study the quantum theory of an 0{N) scalar field on de Sitter geometry at leading order in a 
nonpertative 1/A^-expansion. This resums the infinite series of so-called superdaisy loop diagrams. 
We obtain the de Sitter symmetric solutions of the corresponding, properly renormalized, dynamical 
field equations and compute the complete effective potential. Because of its self-interactions, the 
field acquires a strictly positive square mass which screens potential infrared divergences. Moreover, 
strongly enhanced ultralong-wavelength fluctuations prevent the existence of a spontaneously broken 
symmetry state in any dimension. 
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The impressive observational success of inflationary 
cosmology [H [3] and the quest for high precision measure- 
ments of the cosmic microwave background (CMB) prop- 
erties, such as the power spectrum, or the non-Gaussian 
correlations of temperature anisotropics [3], call for es- 
tablishing a coherent quantum (field) theory (QFT) of 
inflation. Indeed, in this context, the large scale primor- 
dial density fluctuations imprinted on the CMB are of 
purely quantum origin and, as such, are subject to the 
whole complexity of quantum field dynamics. 

The issue of computing QFT radiative (loop) correc- 
tions in curved geometries has a long history [4], which 
has been revived in recent years in the context of infla- 
tionary cosmology [l-Q- In particular, the case of de 
Sitter (dS) geometry has attracted a lot of interest 8- 
[Tbl] . first because of its direct relevance to infiationary 
physics and second, because its high degree of sym- 
metry makes it somewhat simpler [11]. It appears that 
perturbative calculations of radiative corrections of quan- 
tum field dynamics in expanding space-times (ST) share 
many features with similar calculations in flat ST at finite 
temperature ~ due to the presence of an event horizon - 
and out of equilibrium - due to the cosmic expansion. 

For instance, loop diagrams in cosmological STs ex- 
hibit secular terms, which grow unbounded with the 
number of e- folds [g, Q . Those are well-known generic 
features of perturbative calculations for nonequilibrium 
systems [l^ • Another standard problem of perturbation 
theory is the occurrence of infrared (IR) divergences for 
massless fields minimally coupled to gravity, relevant for 
inflation, akin to what happens for light bosonic fields 
at high temperatures or near a second order phase tran- 
sition. Such singularities signal deep deficiencies of the 
approximation scheme and need to be summed [isj . 

Efficient resummation techniques have been developed 
over the years to deal with these issues in equilibrium and 



nonequilibrium QFT in fiat ST. Prominent examples in- 
clude the dynamical renormalization group (DRG) [l3 |. 
or two-particle-irreducible (2PI) functional techniques 
[l5| . First attempts to apply such methods to the 
cosmological case include the work of Riotto and Sloth 
using 2PI and large- techniques [l3l but eventually 
relying on the effective stochastic approach of Ref. [l^ , 
or the DRG analysis of Ref. [l3| • One of the main results 
of these works is the fact that, as was discovered earlier 
by Starobinsky and Yokoyama in the stochastic approach 
p^ , the self-interacting scalar field acquires a curvature- 
induced effective mass, which prevents possible IR diver- 
gences. This is analog to the temperature-induced mass 
in thermal QFT [l5|. 

Here, we study an 0{N) scalar field theory on dS in 
the large- iV limit from the basic principles of QFT. This 
resums the infinite series of so-called superdaisy (or cac- 
tus) loop diagrams and is known to resum high temper- 
ature IR divergences as well as to successfully describe 
the phase structure of 0{N) scalar theories in flat ST 
p^ . We solve the theory in the dS symmetric state con- 
sistent with basic renormalization requirements, the so- 
called Bunch-Davies vacuum, study the phase structure 
of the theory and compute the effective potential in 3 -I- 1 
dimensions [13]. A more extended treatment, including 
the case of a discrete Z2 symmetry (TV = 1) will presented 
elsewhere [2]| . 

We consider a 0(A^)-symmetric scalar field theory with 
classical action (a sum over a = 1, . . . , iV is implied) 



S[v] = - 



' (1) 

with the invariant measure = J d'^'^^x ^/—g, on a d+1- 
dimensional expanding dS background with flat spatial 
sections. In terms of comoving spatial coordinates x and 
conformal time — 00 < 77 < 0, the line element reads 



ds ~ a (rj) idrj — dx ■ dx 



(2) 
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where 0(77) = 
H^Tj^{d^-{d 



-l/Hrj. Here, U = d{d + l)H^ and □ 
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In the large- limit, the field equation reads [13, 

[a + M^{x)]ipa{x)^0 (3) 
where the self-consistent square mass is given by [2^ 



A 



A 



M\x) = + d{d + 1)^H^ + -r (x) + x), (4) 
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in terms of the one-point function {ipa(x)) ~ \/N(t)a{x) 
and the symmetric connected - so-called statistical 
- two-point correlator Fab{x,x') — ^(^(pa{x)(pb{x') + 
^b{x')ipa{x)) with F[x,x') = Faa{x,x')/N. The brack- 
ets represent an average in an appropriate quantum state. 

We wish to investigate dS-invariant solutions of the 
field equations, for which the one-point function (pa is a 
constant and the two-point function F(x, x') — F(s) only 



depends on the dS distance s = (77 +77'' 



2)727777'. 



In particular, F{x, x) = F{1), which implies that the self- 
consistent mass dl]) is constant. 

The statistical two-point function can be written in 
terms of an infinite set of complex mode functions /fc(77) 
in comoving momentum space: 

^^^'^'^ = / -0jd^'''^''^'^'^^iM^)rkiv')} ■ (5) 

Using Eq. ([3|), one finds that the rescaled functions 
Uk{il) — 0(7;)^^ /fc (7;), satisfy the differential equation 

y2_l/4^ 



Uk = 0, 



(6) 



where v = ^Jd? /A — M^/H'^. dS-invariant solutions read 
[23I I: Ufc(7/) = cosh a Sk{ri) + sinh a e''^S'^(r/), where 



Skiv) 



77776 



K'^+3)i/(l)(„fc77), 



(7) 



with Hi^\z) the Hankel function, and where the real 
coefficients a and /3 are k independent. 

The requirement of renormalizability further con- 
strains the possible states of the system. The two-point 
function in the coincidence limit F{x, x) is divergent and 
Eq. Q needs to be regularized and renormalized. In 
principle a dS-invariant regularization (such as e.g. di- 
mensional regularization) should be used. Here, it is suf- 
ficient to use a simple ultraviolet (UV) cutoff on physical 
momenta: k/a{ri) — —kHrj < Ap. Demanding that the 
short-distance (UV) properties of the theory be indepen- 
dent of the ST geometry, one is left with the Bunch- 
Davies vacuum, a = 0, as a unique possibility [23l - [25{ . 

The resummed square mass reads (we make the field 
dependence explicit from now on) 

M^{^^) =m^ + d{d + 1)CH^ + + ^^(^^ (8) 
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with, writing P — Im^v), 

(9) 



We see that this is indeed constant if one uses a physical 
UV cutoff Ap = HA [l^l- A comoving cutoff does not 
lead to a dS-invariant solution. We shall see below that 
this is also crucial to recovering the correct Minkowski 
limit as well as for a proper renormalization procedure, 
involving time- independent bare parameters in Eq. ([I}. 

Now the key observation of the present Letter. It fol- 
lows from the z — behavior: 



-2Rc(i/) 



(10) 

that the loop integral © is only defined if Re(j^) < d/2, 
i.e. the only possible solutions to Eq. (jH) are such that 
M^((/)^) > 0. We emphasize that this is true in any di- 
mension and only follows from the requirement of dS in- 
variance - it is easily checked to hold for any a-vacuum. 
Thus we see that the case of a massless minimally cou- 
pled field, = 0, is not allowed: Whatever the value of 
the tree-level mass, a strictly positive curvature-induced 
mass squared is generated by the field self-interactions. 
We conclude that, in the large- limit, there exits no 
dS-invariant state for an interacting massless minimally 
coupled field 

Let us now examine the consequences of this result. 
Observe first that the field equation 



(11) 



has for unique solution (pa — 0. The strongly enhanced dS 
quantum fluctuations in the far IR, see (1101) , prevent the 
possible existence of a spontaneously broken phase in any 
dimension. This is similar to what happens in fiat ST in 
low dimensions, d < 2 [l^. As we shall see shortly, this 
further implies that any initial symmetry breaking po- 
tential is strongly renormalized by ultralong-wavelength 
fiuctuations, eventually leading to radiative symmetry 
restoration. To discuss these issues in more detail, we 
shall compute the effective potential Vcs{(j)). Recalling 
that dVcs/dcpa = M^i<j)'^)(t>a and exploiting the 0{N) 
symmetry, we have, up to an irrelevant constant Vq: 



V,si<f>) = Vo 



1 



dvM^{v). 



(12) 



Using similar arguments as employed in |26| . it can be 
shown that Vcff{(p) is convex everywhere. 

Before we proceed to actually solving Eq. ([8]), we check 
that we recover the correct flat ST limit _ff — )■ 0, to be 
taken at flxed Ap ~ AH. Note that this also gives the 
relevant equation for the large mass case M'^/H'^ ^ 1, 
for which v = iv k, iM/H. Using the z <C and 



z ^ v behaviors |iJip(z) 



and \H,c{z) 



2e'^'^ /ttVz^+i^ respectively, one recovers the large- iV 
mass equation in the Minkowski vacuum 



M2(02 



A 2 A /■ d'^p 0{Ap~p) 
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(277)^^ 2V?Tm2(02) 



, (13) 



where the contribution from super-Hubble modes is sup- 
pressed by a Boltzmann-like factor ^-^^1"^ with T ~ 
H 1211 the temperature associated with the dS horizon. 
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We now discuss renormalization and specify to d = 3 
spatial dimensions, where the theory is perturbatively 
renormalizable [s^. Using the large-|z| asymptotic ex- 
pansion \Hi^\z)\'' ~ -I [l + (z.2-l/4)/2z2 + 0(z-4)], 
one gets the UV-behavior of the loop integral 



1p 



lnA+21nA+finite (14) 



where "finite" refers to terms which converge in the 
limit A — )> cxD. The first two terms on the right-hand 
side are just the usual Minkowski UV divergences, see 
e.g. Eq. whereas the third term is specific to the 

curved ST. The former are eliminated by standard cou- 
pling and mass renormalization [ltj |, e.g. by introduc- 
ing the renormalized mass and coupling at the arbitrary 
physical scale /x: 1/Ai?(^) — 1/A -I- ln(Ap/^)/487r2 and 
mj^{fj,) / Xr{p) — m'^/X + Ap/487r2. The remaining log- 
arithmic divergence can be eliminated by means of [22| 
12^R{fi)/XR{fi) = 12^/A-hln(Ap///)/247r2. In the follow- 
ing we write the renormalized parameters at the scale 
fi = H as mR{H) = tur and so on. One has the renor- 
malization group (RG) running: 



1 



Aii(M) 



1 



ln(g//i) 
487r2 



(15) 



whereas both 7n\{^) / Xr[^jl) and {^r{p) — \) /Xr{^) are 
RG invariants. Equivalently, (fij^ip) — 2H'^)/Xr{p), 



(fi) + l2^R{fj,)H^ , is RG invariant. 



where fJ-R{fJ.) = itir 

It is interesting to note that the conformal case v = 
1/2, or = 2H^, exhibits an additional symmetry 
which allows one to scale away the cosmic expansion, 
see Eq. (I6|), and which protects the mass from radiative 
corrections as is easily checked using |ffjj2('2^)P = 2/7rz. 
For instance, at = 0, one gets M'^{0) = = 2H'^. 

In the limit M'^/H'^ < 1, Eq. ([8|) can be solved analyti- 
cally. For this purpose, we introduce the small parameter 
e = 3/2 — AP /3H^ ^ 1 and separate the momentum 

integration in Eq. ^ in three parts: = + + 
with k' <^ 1 <^ K <^ Ap. The high-momentum part is 
evaluated by using the large- |z| asymptotic expansion of 
|7?i^^(z)|2. We get, up to the UV divergent terms exhib- 
ited in Eq. (IT4l) . 



rdzz^\Hi'\z)\' ^diY.-- 



1 



To evaluate the intermediate-momentum contribution we 
simply set [H ^ 3/2 and use \hI^/1{z)\'^ = ^{l + z^): 



rdzz'\Hi'H^)\'=- 

Jk' 



K K 

In — + — 



0(e,K 



.l2\ 



(17) 



Finally, the IR contribution is obtained from the z 
expansion \h'^\z)\'^ - {2''T(v) /irz^f [1 + 0(^2)]: 



rdzz'\Hi'\z)f^- 
Jo 



(18) 



where C3 = 2 — 7 — In 2, with 7 the Euler constant. In 
terms of renormalized parameters at the scale jj, — H 
introduced before, we thus get the following equation: 



M\4>^)^^?R■ 



6 ' 



XrH 
487r2 



2 r 



3i72 



M2(02 



2c3 



(19) 



where C3 = C3 -f 1/6 ~ 0.9. The positive solution is [3l| 
M2((/,2) = a(</.2) + Va2(</>2)-t-6 (20) 



with 



487r2 



l2 ' ^ ^ 16^ 



(21) 



The nonanalytic Ai^-dependence is typical of IR effects. 

In this approximation, the effective potential (1121) can 
be analytically integrated and reads (we choose Vb = 0): 

V^s{^) = J- (m4(02) _ ^^2(0)) + M^in^l^5?^. 
2Xr ^ ' ^ 167r2 M2(0) 

(22) 

Full numerical solutions of Eqs. ([5]) and reveal 
that, for a wide range of couplings, Eqs. (I^ni) - (|22l) are in 
fact quantitatively accurate even for nonsmall values of 
AP / . This is because when the infrared term /A/'P 
on the right-hand side of Eq. ()19p is not numerically im- 
portant, the loop corrections are typically suppressed by 
a factor Afl;/487r2. Below, we use Eq. ((20|) to illustrate 
our findings, all of which we checked numerically. 

The above result shows the phenomenon of curvature- 
induced mass generation [l^. Indeed, consider the case 
/i^ -f Ai^(/)2/6 — '32*1, which includes the tree-level mass- 
less minimally coupled field at (/)2 = 0, relevant for infla- 
tion physics, or the Goldstone mode of a tree-level poten- 
tial exhibiting spontaneous symmetry breaking, i.e. with 
/i^ < 0. In this case, the perturbative expansion is ill 
defined due to IR divergences. Self-interactions screen 
these IR effects and generate a nonvanishing effective 
mass and/or nonminimal coupling: 



AP_ 

1p 




(23) 



For weak coupling, where the approximation (jl9l) makes 

,3/2x 



sense, we get AI^/H^ ^ ^/X~^l|A^l-C3XR/^S^l^ + 0{XT2j 



The first term is the Starobinsky-Yokoyama result 18 1, 
obtained here directly from the basic QFT, see also 13 1. 
As discussed in 0, [3 this curvature generated mass 
leads to an enhancement factor H^/AI^ oc I/^/Xr of 
IR-sensitive observables, such as parametrically ©(V Xr) 
non-Gaussianities or corrections to the power spectrum. 

Equation (PO]) also illustrates the phenomenon of ra- 
diative symmetry restoration for /i^ < 0. In that case. 



one gets, aX 4> = 



AP/H^ - A«ir2/167r2|^2^| 



+ 0{Xl). 



Notice the parametric Xr dependence which may lead to 
an even stronger enhancement of IR-sensitive observables 
than mentioned above. As an illustration. Fig. [T] shows 
an exact numerical solution of Eqs. ([S]) and (IT^ . 
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FIG. 1: The effective (red) vs tree-level (blue) potentials as a 
function of for = — 2 and Aij = 10, all in units of H. The 
inset shows the square mass M^(<^^) (red) and its tree-level 
expression ^\ + \r4>^ /6 (blue). Here, the numerical result 
agrees with Eqs. (|20[) - (I22I ) over the whole range of </). 

In conclusion, we find that the strongly enhanced 
ultralong-wavelength scalar field fluctuations in dS pre- 
vent IR divergences by generating a strictly positive 
curvature-induced square mass and, as a consequence, 
forbid the possible existence of a spontaneously broken 
phase. This is a genuine nonperturbative effect. It is an 
interesting question to investigate whether this result re- 



mains true beyond the large-N limit 2l|. Together with 
Refs. [13, this work provides a step toward the pre- 
cise understanding of such nonperturbative effects in dS 
from first principles. 

We emphasize that the phenomenological relevance of 
these results requires a rather long period of (nearly) dS 
expansion in order for such IR fluctuations to develop. In 
the small mass limit ^T9\\ . the relevant wavelength are as 
large as exp(37J^/2M^). This is certainly not a problem 
for e.g. eternal infiation scenarios. 

The phenomena of mass generation or symmetry 
restoration are similar to finite temperature effects in 
fiat ST. It is a very interesting question whether other 
important thermal effects such as damping of unequal 
time correlators [l2| occur in dS. This requires calcula- 
tions beyond the present approximation [111, [l^ . 

Various techniques have been developed to deal with 
IR and secular issues in equilibrium and nonequilibrium 
QFT in flat ST [H H 111. It remains to be inves- 
tigated whether these can help resolving similar ques- 
tions in curved geometries, or if specific methods will be 
needed. 

We acknowledge useful discussions with R. Parentani 
and J. Pawlowski. 
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